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We present a cooling method controlled by a physical cooling radius that defines a scale below which fluctuations 
are smoothed out while leaving physics unchanged at all larger scales. This method can be generally used as a 
gauge invariant low pass filter to extract the physics from noisy MC configurations. Here we apply this method 
to study topological properties of lattice gauge theories where it allows to retain instanton-anti-instanton pairs. 



Since MC configurations are noisy at small dis- 
tances, smoothing is needed to observe the physi- 
cal structm-e. One approach is "cooling" a lo- 
cal minimization of a given lattice action. Cooling 
works as a diffusion process, smoothing out in- 
creasingly large regions ||] . Thereby the physical 
spectrum is affected, in particular the string ten- 
sion drops rapidly with cooling. Topological exci- 
tations may be distorted by bad scaling properties 
of the action and by instanton - anti-instanton 
(I-A) annihilation (pairs are not minima of the 
action). We thus need: 

1) to use an action with practically scale invariant 
instanton solutions and a dislocation threshold to 
eliminate UV noise (dislocations), and 

2) to control cooling by means of a physical scale 
such that no monitoring or engineering which 
may introduce uncertainties is necessary. 
Restricted Improved Cooling (RIG) fulfills these 
requirements |^ . RIG preserves physics at scales 
above a cooling radius r which can be fixed un- 
equivocally beforehand, while smoothing out the 
structure below r. In particular the string ten- 
sion is preserved beyond r, instantons are stable, 
dislocations are eliminated, and I-A pairs are re- 
tained above a threshold defined by r. 

Properties of RIC. RIG uses the action of the 
Improved Cooling algorithm (IG) with 5 planar, 
fundamental Wilson loops m . This action is cor- 



rect to order 0{a^) and has a dislocation thresh- 
old pq ~ 2.3a, below which short range topologi- 
cal structure is smoothed out (note that pQ ^ 
in approaching continuum). Above po, instan- 
tons are stable to any degree of cooling (how- 
ever, I-A pairs annihilate). The corresponding 
improved charge density using the same combina- 
tion of loops leads to an integer charge already 
after a few cooling sweeps and stable thereafter. 

Here we shall restrict ourselves to SU{2), for 
the general case see Recall that the cooling 
algorithm is derived from the equations of motion 
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where W is the sum of staples connected to the 
link U^{x) in the action. It amounts to substitute 
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|2 ^ ^Tr{WW^). (2) 



* Partial support from the European Networfc "Finite Tem- 
perature Phase Transitions in Particle Physics" . 



We define RIC by the constraint^ that only those 
links be updated, which violate the equation of 
motion by more than some chosen threshold: 

C/ -> y iff A^{xf = ^Tr(l - UV^) > S'\ (3) 

We have Al{x) oc -Tr((D,,F^^(a;))2) in contin- 
uum limit (5|. Thus 6 controls the energy of the 
fluctuations around classical solutions and acts 
as a filter for short wavelengths. Since it uses the 
same action RIG has the same scaling properties 

^We thank F. Niedermayer for this suggestion. 
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as IC. However, since RIC does not update links 
already close to a solution, it changes fewer links 
after every iteration until the algorithm saturates. 

defines a constrained minimization and the 
smoothing is homogeneous over the lattice. 

The next step is to relate the parameter S which 
defines the cooling to a physical scale. For Yang 
Mills theory the latter should involve the string 
tension a. We calculate the "effective mass" M{t) 
from correlation functions of spatial Polyakov 
loops separated by t steps in time. Asymptoti- 
cally M{t) Nga up to finite size corrections. 

For the calculations we generate configurations 
using the Wilson action on the lattices: 

12^ X 36, p.b.c, 13 = 2.4 (a = 0.12 fm): 800 
configurations with 100 sweeps separation; 

24"*, twist in time, /3 = 2.6 (a = 0.06 fm): 350 
configurations with 200 sweeps separation. 
We do 20000 thermalization sweeps. We use S = 
2.89, 4.92, 11.57, 23.15 and 46.30 fm~3. 

For illustration we present in Fig. |l| M{t) for 
the 24^ lattice. Defining r{d) as the distance t 
at which M{t) on 5-cooled configurations starts 
to agree with the uncooled value (obtained by 
fuzzing and fitted to a smooth function of t) we 
arrive at the results in Fig. ^, compatible with 

r{S) ~ 0.8 S-^/^, (4) 

showing the correct scaling behaviour, as ex- 
pected since S has a continuum limit. 
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Figure 1. M(t) for the a = 0.12 fm lattice. The 
horizontal band represents standard results for a. 

The behaviour of instantons is similar under 
RIC and IC whereby the former has a somewhat 
smaller, (5-dependent dislocation threshold po('^)- 
I-A pairs, however, do not annihilate under RIC 
the way they do for the other cooling methods 
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Figure 2. The smoothing scale of RIC: r(S) 
vs (5^^/"^ for the a ~ 0.12 fm (diamonds) and 
a — 0.06 fm (squares) lattices (the middle points 
are slightly displaced around the true positions for 
readability). Also shown is the fit 

(including IC): since the distortion of the part- 
ners in a pair depends on the overlap, RIC pre- 
serves pairs below some overlap-threshold. More 
precisely, it turns out that there is a well defined 
relation between A and S}^ = IGn^ — S^^, there- 
fore RIC stabilizes I-A pairs with 5int below a 
threshold which is a function of S, hence of r. 

SU(2) topological properties by RIC. The 

total, improved topological charge stabilizes to 
integer values after only few IC or RIC sweeps. 
Therefore the topological susceptibility appears 
very stable if S is small enough to eliminate dis- 
locations, shows correct scaling and agrees with 
the Witten-Veneziano relation. 

To describe the details of the instanton ensem- 
ble one needs to recognize the particular I's and 
A's. Since the objects can be very distorted, our 
description relies on two assumptions: a) instan- 
tons should appear as local self-dual peaks in the 
action and charge density, and b) only pairs with 
Sj^ considerably smaller than 167r^ should be 
considered as such. We approximate the action 
and charge density by a superposition of I's and 
A's parameterized through the BPST formula for 
an instanton of size p. We measure the adequacy 
of this ansatz through the fraction of the total 
action (charge) reproduced in this way. We rep- 
resent the data in physical units using the corre- 
sponding a and the r{S) obtained from the string 
tension analysis without any further fit or tuning. 

The main results of our analysis are: 
1. The density of instantons increases drastically 
with decreasing smoothing scale r. 
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Figure 3. Normalized size distributions on the 
a = 0.12 fm lattice (squares, triangles) and a — 
0.06 fm lattice (crosses). Vertical dotted (dashed) 
lines indicate the IC dislocation threshold po = 
2.3a for a = 0.12 (0.06) fm. 




Figure 4. Non-normalized overlap distributions 
for a = 0.06 fm and 6 = 2.89 (circles), 4.92 
(bars), 11.57 (squares) and 23.15 (crosses) fm~^ . 

2. Size distributions depend only weakly on r. 

3. The typical I-A distance diA (from an I to the 
nearest A) seems given by the smoothing scale. 

4. The overlap {pi + pA)/2djA increases strongly 
with decreasing r and the fit quality deteriorates. 
All these results are largely invariant under 
rescaling of the lattice spacing a. They are sum- 
marized in the Figures |3[ || and ||. 

It appears therefore that there are topological 
features of the vacuum structure, like suscepti- 
bility and charge distributions and to a fair de- 
gree also size distributions, which are well defined, 
show correct scaling and are rather independent 
on the smoothing scale once the noise has been 
damped sufficiently. However, I(A)-density, over- 
lap and I-A distance distributions are not of this 
kind: although these quantities seem to scale cor- 
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Figure 5. Dependence of the average size (tri- 
angles), the average I-A distance (diamonds) 
and the instanton density (squares) on the cool- 
ing radius r. The dashed curve is {N/V) = 
(A^/V^)asymptotic + C / r"^ . Open (filled) symbols 
correspond to the a = 0.12 (0.06) fm lattice. 

rectly with the cut off, their behaviour with r sug- 
gests that looking at smaller and smaller scales a 
continuous spectrum of fluctuations emerges and 
a description in terms of I's, A's and I-A pairs 
becomes less and less meaningful. This may ex- 
plain both the agreements and the disagreements 
between various analyses - see, e.g. 

Generally RIG proves itself a good instrument 
for probing different scales in a controlled way 
and extract physics from noisy configurations. 
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